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Abstract 

A generalized continuum framework for the theory of solute transport in fluids is 
proposed and systematically developed. This framework rests on the introduction 
of a generic force balance for the solute, a balance distinct from the macroscopic 
momentum balance associated with the mixture. Special forms of such a force 
balance have been proposed and used going back at least as far as Nernst's 1888 
theory of diffusion. Under certain circumstances, this force balance yields a Fickian 
constitutive relation for the diffusive solute flux, and, in conjunction with the solute 
mass balance, provides a generalized Smoluchowski equation for the mass fraction. 
Our format furnishes a systematic procedure for generalizing convection-diffusion 
models of solute transport, allowing for constitutive nonlinearities, external forces 
acting on the diffusing constituents, and coupling between convection and diffusion. 

1 Introduction 

The transport of a solute in an incompressible fluid mixture of density g is commonly 
modeled by the convection-diffusion equation^] 

c = div (D(c) grade) , (1-1) 

with c the solute mass fraction and D > the mass fraction-dependent diffusivity. This 
equation is often derived (see, for example, Landau & Lifshitz de Groot & Mazur [Q, 
and Cussler [||) by adjoining to the solute mass balance 

qc = -div j, (1.2) 

a constitutive relation 

3= -QD{c) grade (1.3) 

of the sort proposed by Fick B for the diffusive solute flux j. 

Here, we propose and systematically develop an alternative approach to solute trans- 
port. Our approach is based on the viewpoint that to each independent kinematic process 
there should correspond a system of power- conjugate forces subject to a distinct force 



1 Here, the operations of material and spatial time-differentiation are denoted with a superposed dot 
and a superposed grave accent, respectively; thus, for instance, c = c+v-grad c, with v being the mixture 
velocity. 
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balance.^ Specifically, in the context of solute transport, we identify mixture convec- 
tion and solute diffusion as independent kinematic processes: mixture convection is a 
macroscopic phenomenon described kinematically by the convective mixture velocity v, 
whereas solute diffusion is a microscopic phenomenon described kinematically at a con- 
tinuum level by the diffusive solute flux j, an object that measures the motion of solute 
particles relative to the macroscopic motion of the mixture comprised of solute and sol- 
vent. Further, to account for power expenditures associated with the diffusive motion 
of atoms, we introduce at a continuum level a system of forces conjugate to the diffu- 
sive solute flux j. Since they are associated with the microscopic processes governing 
the diffusion of the solute, we use the expression "microforce" to identify these forces. 
Additionally, we require that these microforces be consistent with a solute microforce 
balance. This balance is a postulate ancillary to the conventional momentum balance, 
which involves the forces that act power conjugate to the convective mixture velocity v. 
Importantly, we do not assume that the diffusive solute flux j is determined by a con- 
stitutive relation — Fickian or otherwise. Rather, we treat j as an independent kinematic 
variable whose evolution is governed by the microforce balance. 

There is ample historical precedence for basing solute transport on additional, inde- 
pendent force balances. Nernst following upon van't Hoff's [{[7J analogy between the 
osmotic pressure w of the molecules of a sufficiently dilute solute and the pressure of an 
ideal gas, viewed diffusion as driven by osmotic pressure gradients and developed a theory 
of diffusion based on a balance of forces acting on the solute molecules. Nernst proposed 
such force balances to account for the influences of osmotic, viscous, electrostatic, gravi- 
tational, centrifugal, and magnetic forces. For example, for electrolyte solutions, Nernst 
assumed separate force balances for the positive and negative ions. Subsequently Planck 
pl| also exploited these ideas in his work on electrolyte solutions. Further, in develop- 
ing his theory of Brownian motion, Einstein |l9|, ||l|] employed force balances on the 
Brownian particles. Neglecting the inertia of the diffusing solute, we may write Nernst's 
solute force balance as 

1 M 

h + £ = — grad w — — grad w, (1.4) 
n qc 

wherein h is the hydrodynamic drag force acting on the solute particles during motion 
through the solvent, £ represents the external forces acting on the solute particles, M 
is the molecular weight of the solute, and n is the number density — so that the right- 



hand side of (1.4) can be interpreted as the osmotic force per solute particle. Typically 
the hydrodynamic drag force is taken to be linear in the solute velocity relative to the 
mixture, whereby 

h = -£j, (1.5) 

with £ a non-negative coefficient called the frictional resistance and C _1 called the mo- 
bility^ For an ideal solution, the osmotic pressure w is linear in the number density n, 
that is 

gk B 9c 

w = k-Qvn = — — — , (1-D) 
M y ' 
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A foundation for this viewpoir 



t can be devised using the principle of virtual power, see, for example, 



Germain g], Antman & Osborn [ jj, Maugin and Fremond 0|. We follow here an alternative, but 
essentially e quiv alent, approach taken by Erickscn M, Goodman & Cowin [nXf, Capriz & Podio-Guidugli 
jn|, Capriz [0, Fried & Gurtin |Q, Fried and Gurtin [[lij. 

3 Since we take the solute flux j as a generalized velocity, our definitions of frictional resistance and 
mobility differ from the usual ones by a factor of qc. 
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with &b being Boltzmann's constant and 9 being the absolute temperature. Thus, in an 
isothermal, constant density mixture, 

gradw — k^8 gradn — ^—j^j— grade, (1-7) 



which, in conjunction with (1.4) leads to the flux 



k B , £ , 
3 = -^-gradc+-. 1.8 



Hence, prov ided t hat external forces are absent, the balance (1.4) and constitutive re- 
lations ( |l.5| ) and (1.7) yield an expression similar to the Fickian constitutive equation 
( |l.3| ) with the classical relation (see Nernst Page 615]) 

between diffusivity and frictional resistance (or mobility). 

There is, however, disagreement over the form and interpretation of the solute force 
balance. For example, Einstein |^2|, § 4] asserted that the osmotic pressure is a "virtual" 
force, so that it cannot be treated as a force acting on the individual molecules^] In fact, 
Einstein |]20f , in his original treatment of Brownian motion, proposed instead a fictitious 
force k that acted on the individual particles but that balanced the osmotic pressure 
gradient: 

M 

k = — gradw. (1-10) 

QC 

Einstein j2(J] then took the solute force balance as 

h + l = k (1.11) 

and thereby claimed (see Einstein |22]) to have resolved the apparent difficulty of having 
a virtual force act on the individual particles. 

Additionally, Debye |2^, Page 106] asserted that osmotic forces are not "true" me- 
chanical forces, so that the solute force balance should be instead 

h + £ = 0. (1.12) 

In this alternative point of view, the flux j is decomposed as a sum 

3 = 3d+3f (1-13) 

involving a purely diffusive part j D (usually called the "ordinary diffusion" ) and a part 
j F due to external forces (usually called the "forced diffusion"). Commonly, Fick's first 
law is invoked for the ordinary diffusive flux, so that 

j D = -Dgradc, (1.14) 

and the flux due to the external forces is assumed to have the form 

3f = ~ c (1-15) 



4 The exact meaning of this statement is unclear to the authors. 
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which follows from (1.5) and ( 1 . 12f ) . Combining ( 1.15 ) with the balance (1.2) leads to 



the generalized convection-diffusion equation 

c = div (D grade -), (1-16) 

sC 

typically referred to as a Smoluchowski |24| equation. The diffusivity D can be deter- 
mined by considering the special case where the external forces are given by a potential, 
viz., 

£=-gradV, (1.17) 

and requiring the steady-state solution to correspond to that of a Boltzmann J2q| distri- 
bution, so that 

/-V\ 

c = c exp(— ). (1.18) 



This yields the relation (1.9), which is usually enforced even when ( [1.18 ) no longer holds. 



This latter format was used by Debye |26|, |27| in his theory of rotary diffusion in 
dielectrics to include the effect of the external electric field on the diffusing dipolcs. 
It is still commonly used to obtain extensions of the Fickian constitutive relation to 
include, for instance, the effects of external electric fields on electrolytes. Additionally 
it is used to model the interaction of diffusion and external fields on dilute suspensions 
(see, for example, Giesekus [p8[ and, more recently, Doi & Edwards p9fl). However, Bird, 



Armstrong, & Hassager [|30[ based their approach to polymeric suspensions on a force 
balance such as ( |1.4| ) that explicitly includes a contribution attributed to the Brownian 
forces on the polymer molecule (see, for example, |R)L Equation 13.2-1]). 

Interestingly, in the case of linear constitutive relations, the approaches described 



above lead to the same resulting equations — namely (1-16) and (|1.9|). However, the 



assumptions involved are rather different. In particular, it is not clear how to general- 



ize either approach to allow for constitutive nonlinearities. For example, in ( 1.13 ) the 
ordinary diffusive flux j D is determined by a constitutive relation, whereas the forced 
diffusive part j F is determined by a balance law — a procedure highly atypical for a con- 



tinuum theory. Furthermore, the validity of (1.9), derived from an assumed steady-state 
solution, is not clear in unsteady situations. 

In our opinion, these issues can be resolved by a more general treatment of the 
notion of forces in a binary mixture. Consequently, in the following, we develop a the- 
oretical framework for diffusion based on a generic microforce balance without a priori 
assumptions on the nature of the microforces or of steady-state concentration profiles. 
Furthermore, we allow for constitutive nonlinearities and external fields. We follow the 
procedure of modern thermodynamics, wherein consistency with the second law restricts 
the form of the constitutive relations. Our proposed format provides a systematic proce- 
dure for the generalization of convection-diffusion equation models for solute transport. 
In particular, it leads to a force balance that explicitly includes the diffusive or Brownian 
forces acting on the solute. Only in special cases, such as linear constitutive relations, 



can this force balance be expressed in the oft used form (1.12) with the corresponding 



decomposition (1.13) of the flux proposed by Debye. 

For clarity, we restrict our attention to purely mechanical processes in a binary mix- 
ture of solute and solvent, ignoring thermal and other effects. In this context, the first 
and second laws of thermodynamics are replaced by a free- energy imbalance that ac- 
counts for power expenditures associated not only with mixture convection but also with 
solute diffusion. That diffusion does, indeed, lead to the expenditure of power was clearly 
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recognized and utilized by van't HofT |17| as an experimental means of measuring the 
osmotic pressure. First, we consider simple diffusion without mixture convection. This 
special case allows us to develop the theory without the complications of coupling with 
the macroscopic mixture motion. We then treat both convection and diffusion. 



2 Simple theory accounting only for diffusion 

We consider a binary mixture of solute and solvent that occupies a fixed region 1Z of three- 
dimensional space. In the present context, where we ignore the macroscopic motion of 
the mixture, the mixture density is constant. Without loss of generality, we take that 
constant to be unity. 

2.1 Basic laws 

The simple theory is based upon the following basic laws: 

• solute mass balance; 

• solute microforce balance; 

• mixture free-energy imbalance. 

We formulate these laws in global form for an arbitrary subregion V of TZ. The outward 
unit-normal field on the boundary dV of V is denoted by v. 

2.1.1 Solute mass balance 
We introduce fields^] 

c solute mass fraction, 
J diffusive solute flux, 
m external solute supply, 

in which case the integrals 

cdv, J jv da, and J m dv (2.1) 

v dv v 

represent, respectively, the solute mass in V , the solute mass added to V , per unit time, 
by diffusion across dV, and the solute mass added to V , per unit time, by external 
agencies. 

Solute mass balance is the postulate that, for each subregion V of TZ and each instant, 



cdv=— J j-vda + J m dv. (2-2) 

V dV V 

Equivalently, since V is arbitrary, we may enforce solute mass balance via the local field 
equation 

c=—divj + m. (2-3) 



5 Since the mixture density is assumed to be unity, there is no distinction between solute mass and 
solute mass fraction, or between solute mass flux and solute mass fraction flux. 
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2.1.2 Solute microforce balance 

Associated with the evolution of the diffusive solute flux j, we introduce a system of 
power-conjugate microforces, consisting of 

II solute microstress tensor, 

h internal solute body microforce, 

£ external solute body microforce, 

in which case the integrals 

III/ da, J hdv, and J i dv (2.4) 

TV V 

represent, respectively, the solute microforces exerted on the region V by the solute 
microtraction distributed over dV, by agencies within V, and by agencies external to V , 
respectively. In the present mixture context, internal forces arise from the action of the 
solvent on the solute, or from solute-solute interactions. External forces are those arising 
from sources remote from the mixture, such as those arising from applied electromagnetic 
or gravitational fields. The microstress tensor 77 allows for surface tractions that are not 
parallel to the surface normal. We assume that the inertia associated with solute diffusion 
is negligible and, hence, impose a solute microforce balance, which is the postulate that, 
for each subregion V of 1Z and each instant, the resultant microforce acting on V must 
vanish: 

J nvda + J(h + £) dv= 0. (2.5) 

av v 

Equivalently, since V is arbitrary, we may enforce solute microforce balance via the local 
field equation 

divll + h + l = 0. (2.6) 



A comparison to the Nernst force balance (1.4) suggests an interpretation of the quantity 



U as a generalized stress tensor that allows for the possibility of tangential osmotic forces 
on a semipermeable membrane. At this point, we impose no constitutive assumptions 
concerning the microforces. In particular, we allow for the possibility that II vanish 



identically, which would correspond to the alternative balance ( |l,12[) proposed by Debye. 



2.1.3 Free-energy imbalance 

In this purely mechanical theory, the first and second laws of thermodynamics are re- 
placed by a free-energy imbalance. To formulate this imbalance, we introduce fields 

ip mixture free-energy density, 
fi solute diffusion potential, 

in which case that the integrals 

ijj dv and J /im dv (2-7) 
v v 

represent, respectively, the free-energy in V and the rate at which free-energy is added 
to V by the external supply of solute to V . 
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Further, the integrals 

Ilujda and / i-jdv (2. 



av v 

provide an accounting of the power expended on V by the solute microtraction distributed 
over dV and by the agencies acting external to that region. These expressions represent 
the power expended by the diffusion of the solute particles. 

Consistent with our assumption that solute inertia is negligible, we neglect the kinetic 
energy of solute diffusion. Free-energy imbalance is, then, the postulate that, for each 
subregion V of 1Z and each instant: 



J ipdv < J nvjda + J + fim) dv. (2.9) 
v av v 

Equivalently, since V is arbitrary, we may enforce free-energy imbalance via the local 
field inequality 

-■tji + fic-h-j+ill + /j,l)-gradj > 0. (2.10) 



2.2 Constitutive equations 

The preceding balance laws do not form a determinate system of equations. We now 
develop constitutive equations that will provide closure relations. We assume that the 
free-energy density tp, the solute diffusion potential /i, the solute microstress tensor 77, 
and the internal solute body microforce h are determined constitutively as functions of 
the solute mass fraction c and the diffusive solute flux j, so that 

(V, M, n, h) = $(c,j), A(c,j), n(c,j), h(c,j)). (2.11) 

We emphasize that, in contrast to the standard approach to the theory of mass diffusion, 
the diffusive flux j is not given by a constitutive relation. 



Inserting the constitutive relations (2.11) into the local free-energy imbalance (2.10), 
we arrive at the functional inequality 

(^(c,j)-A(c ) j))c+^(c,j)-i + Mc,j)-j-(l7-(c,j) + A(c,j)i)-gradj<0, (2.12) 

and, following the procedure founded by Coleman & Noll j3l[] in their incorporation of 
the second law into continuum thcrmomcchanics, conclude that: 

(i) the constitutive response functions tp, fi, and U delivering the free-energy density 
tjj, the diffusion potential fx, and the stress II must be independent of the diffusive 
mass flux j, witbQ 

tp = i)(c), fi = i)'(c), and n = -tp'(c)l; (2.13) 



(ii) the constitutive response function h for the internal solute body microforce h must 
be consistent with the residual inequality 

h(c,j)-3<0. (2.14) 

6 Here, a prime indicates with respect to the argument of a function depending on a scalar field. It 
is noteworthy that the solute microstress reduces to an isotropic tensor because we do not include the 
gradient grad j of the diffusive solute flux j in the list of independent constitutive variables. 
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Granted smoothness of the response function h, a result appearing in Gurtin & 
Voorhees [^|| yields the general solution 



h(c,j) = -Z{c,j)j 



(2.15) 



of (2.14), where the symmetric generalized frictional resistance tensor Z must obey 

j-Z(c,j)j>0. (2.16) 
Objectivity further requires that under the transformation 

ct-^c, j^Qj, (2.17) 
with Q an arbitrary rotation, the quantities ip, fj,, II, h, and i transform as 

ip^ip, fi^fi, II^QIIQ\ h^Qh, £^Q£, (2.18) 
which, in particular, implies that 

Z(c,3) = ((c,j)l, 3 =\j\, (2.19) 
where the scalar frictional resistance C must obey 

C(c,j)>0 (2.20) 
for all (c,j). Combining and ( pJ9| ), we see that 

fr=-C(c,j)j, (2.21) 

so that the internal solute body microforce h must be collinear with the diffusive solute 
flux j. 

The foregoing results show that the behavior of a medium of the sort considered here 
is completely determined by the provision of two constitutive response functions: (j) ip 
determining the mixture free-energy density as a function of the solute mass fraction 
c; and, (ii) the scalar frictional resistance £, which determines the internal solute body 
microforce h and, in general, may depend on both the solute mass fraction c and the 
magnitude j of the diffusive solute flux j. 

The governing equations that arise on substituting the foregoing thcrmodynamically 
consistent constitutive equations in the local field equations (2.3) and (2.6) expressing 
solute mass balance and solute microforce balance are 



c + div j = m, 
C(c,j)j + grad (ft(c))=t. 



(2.22) 



Provided that £ is nonvanishing, (2.22)2 yields an expression 



C(c,j) 



grad c 



C(c,j) 



(2.23) 



for the diffusive solute flux j. This leads, in conjunction with ( 2.22 )i, to the generalized 
Smoluchowski equation 



, /V"(c) , I ■ 



(2.24) 
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Note that, when the frictional resistance £ depe nds on j, the evolution equation ( 2.24 ) 
does not decouple from the microforce balance (2.22)2- In this case, the flux j does not 
separate into distinct ordinary and forced diffusive contributions. Comparing ( 2.23| ) with 
the conventional relation (1.3), we obtain a generalized expression 



D(c,j) 



C(c,j) 



(2.25) 



for the diffusivity, which is not predicated on a particular steady-state solution and, 
moreover, is not restricted to the linear regime. Hence, within our framework, the Fickian 
constitutive relation ([lj]) can be interpreted as an expression of the solute microforce 
balance, granted that the internal body microforce is linear in the diffusive solute flux 
and that the external solute body microforce vanishes. In the special case governed by 



(1.17), the steady-state solution is 



(j, = i/>'(c) = — V + constant, 



(2.26) 



which need not correspond to the Boltzmann distribution ( 1.18 ) 



Our results for the microforce bala nce are consistent with the Nernst balance (1.4). 
In fact, a comparison of (1.4) to ( 2.22^ 2 shows that they are identical in form if 



w'(c) = —4>"(c) = — u'(c). 



(2.27) 



Thus, (2.13) shows that the osmotic traction acts only normal to a semipermeable mem- 
brane and is determined by the mixture free-energy density. Note that (2.22) 2 is consis- 
tent with the Debye ( 1.12 ) balance only when the frictional resistance £ is independent 
of the flux j. 



2.3 Special case: diffusion of non-interacting spheres 

We now specialize the preceding results to the case of diffusing, non-interacting spher- 
ical particles, which corresponds to the dilute Brownian spheres originally treated by 
Einstein fl9|| . Our framework requires only the specification of constitutive relations for 
the mixture free-energy density and the frictional resistance. We stipulate that ip be 
determined by a response function ip of the form 

$(c) =fc B 0elogc, (2.28) 

which corresponds to an ideal mixture in the sense that it consists (to within an arbitrary 
additive constant) only of a classical entropic contribution, a contribution that drives 
diffusion. This choice does not encompass particle-particle interactions. Insertion of 
( |2.28| ) into ( |2.27| ) leads to 



Further, we require that the response function h be linear in the diffusive solute flux j 
with the specific form 

h = h(c, 3 ) = -C(c) 3 , with C (c) = ^, (2.30) 

so that the internal body microforce is consistent with Stokes' [^3| law for the drag on 
a rigid spherical particle slowly translating in an incompressible Newtonian fluid. Here, 
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rj s is the shear viscosity of the pure fluid solvent and r is the particle radius. Finally, 
we suppose that the external body microforce is generated by the gradient of a potential 
according to ( 1.17 ). 

Granted these assumptions, we may use ( [2.23 ) to obtain 



3 = 



■ grad c — grad V, 



Qirrj s r 



6irrj s r 



(2.31) 



which, when inserted into the mass balance ( 2.22 )2 yields the generalized diffusion equa- 
tion 



c = D div (grad < 



; grad V) , 



with the diffusivity D having the explicit form 

k B 



D = 



6irrj s r 



(2.32) 



(2.33) 



consistent with the diffusivity expression originally derived by Einstein fl9|| for Brownian 
spheres. Thus, the contribution of ip in ( [2.31 ) leads to a classical, linear second-order 
diffusion term that tends to randomize concentrations, which justifies our description of 
it as entropic. Additionally the steady-state mass fraction has the classical Boltzmann 
form Q. 



Importantly, our derivation of (2.32) and (2.33) relies neither on a Fickian constitutive 
relation for the diffusive mass flux nor on an assumed form for the mass fraction distri- 



bution in equilibrium. Rather, the relation (2.33) is a consequence of the thermodynamic 
theory and the specialized ideal constitutive expressions (2.28) and (2.30). 



3 Theory accounting for both mass diffusion and con- 
vection 

We now account for motion of the mixture as described by a convective velocity v. 
Further, we allow the region 1Z of three-dimensional space occupied by the mixture to 
evolve with time. 

3.1 Basic laws 

The theory is based upon the following laws: 

• mixture mass balance; 

• solute mass balance; 

• mixture momentum balance; 

• moment of mixture-momentum balance; 

• solute microforce balance; 

• mixture free-energy imbalance. 

As in our treatment of the theory without mixture convection, we formulate the basic 
laws in global form for an arbitrary material subregion V of 1Z. 
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3.1.1 Mixture mass balance 

We now allow for convection and variations in the mixture mass density, so that we 
introduce fields 

q mixture mass density, 
v mixture velocity. 

Mixture mass balance is the postulate that, for each material subregion V of 1Z and 
each instant, 



gdv=0. (3.1) 



v 



Equivalently, since V is arbitrary, we may enforce mixture mass balance via the local 
field equation 

g + gdivv = 0. (3.2) 

3.1.2 Solute mass balance 

With mixture convection taken into account, the global statement of solute mass balance 



for a material subregion V of 1Z has the form (2.2) considered in the absence of convection 



provided that c is replaced by gc. Thus, bearing in mind the mixture mass balance ( p.2[ ), 
the local equivalent of solute mass balance now reads 

gc=—divj + m. (3-3) 
3.1.3 Solute microforce balance 



With mixture convection taken into account, the global and local statements, (2.5) and 



( |2.q ), of solute microforce balance remain unchanged from those introduced earlier. 

3.1.4 Mixture momentum balance and moment of mixture- momentum bal- 
ance 

Associated with the convective mixture velocity, we introduce fields 
T mixture stress, 

b external mixture body force density, 
in which case the integrals 

Tv da, and / b dv (3.4) 



represent the mixture forces exerted on the region V by the mixture traction distributed 
over dV, by agencies within V , and by agencies external to V , respectively. 



In addition to the mixture mass balance (3.2), the solute mass balance (3.3), and the 



solute microforce balance (2.6), we enforce balances of mixture momentum and momen- 
tum of mixture momentum. 
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Mixture momentum balance is the postulate that, for each material subregion V of TZ 
and each instant, the time-rate at which the convective momentum within V changes be 
equal to the resultant of the convective forces acting on V: 



J gv dv= J Tu da + J b dv. (3.5) 



V dV 



Further, moment of mixture momentum balance is the postulate that, for each material 
subregion V of TZ and each instant, the time-rate at which the moment of the convective 
momentum within V changes be equal to the resultant of the convective torques acting 
on V: 



J xAgvdv= J xATvda + J xAbdv. (3.6) 

V dV V 

Equivalent to these balances, we have the local field equations 

qv = divT + b and T = T T . (3.7) 

3.1.5 Free-energy imbalance 

Our previous statement of free-energy imbalance must be modified to account for vari- 
ations in the mixture mass density, the kinetic energy of the mixture, and the power 
expended by the mixture stress and external mixture body force. Specifically, we now 
require that for each material subregion V of TZ and each instant: 



J g(ip + ±M 2 ) dv < J (n T j + T T v)-L> da + J (£-j + fim + b-v) dv. (3.8) 

V dV V 

Equivalcntly, since V is arbitrary, we may enforce mixture energy imbalance via the local 
inequality 

-gip + qhc - h j + T D + (n + [il)-gradj > 0, (3.9) 

with 

D = \(gradv + (gradv) T ) (3.10) 

the symmetric strain-rate. 

3.2 Constitutive equations 

We consider the cases of compressible and incompressible mixtures separately. 
3.2.1 Compressible mixture 

To the lists of dependent and independent constitutive variables we add the convective 
stress T and the mixture density g and strain-rate D, respectively, so that 

(V, V, n, h, T) = (j,(z), fi{z), fl(z), h(z), T(z)), (3.11) 

where, for brevity, we have introduced 

z = {g,c,j,D). (3.12) 
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Inserting the constitutive relations ( 3.11 ) into the local energy imbalance ( |3.9| ), we 
arrive at the functional inequality 

dip dip dip ~ 

e(-fc( z ) ~ K z ))c + Q-g-(z)-j + Qg^(z)-D + h(g, c,j,D)j 

- (T(z) -g 2 ^(z)l)-D- (fl(z) + fi(z)l)-gradj < 0. (3.13) 

Hence, proceding as in the our treatment of the theory without mixture convection, 
we find that: 

(i) the constitutive response functions ip, ji, and 17 delivering the mixture energy 
density ip, the mixture diffusion potential /i, and the solute microstress U must be 
independent of j, with 

d - 

ip = ip(g,c), fi= —i/;(g,c), and JT = -fi(g,c)l; (3.14) 

(ii) the constitutive response function h for the internal solute body force density h 
and the mixture stress T must be consistent with the residual inequality 

h(g,c,3,D)-j-T vis (g,c,3,D)-D < 0, (3.15) 

in which 

T vis (g,c,j,D) = g 2 -^-ijj(g,c)l + T (3.16) 

determines the dissipative contribution T v ; s to the mixture stress. 

Further, granted smoothness of the response functions h and T v j s , the previously 
employed result of Gurtin & Voorhees |52| yields a general solution 



h = -Z(g,c,j,D)j-Y(g,c,j,D)D,' 
T vis = U(g,c,3,D)j + V(g,c,3,D)D^ 



(3.17) 



of (3.15), where Z, Y, U, and V are isotropic due to objectivity and must obey 

3-Z(g, c,3, D) 3 + 3 -Y{g, c, 3 , D)D + D-U(g, c, 3, D) 3 + D-V(g, c, 3, D)D > (3.18) 

for all (g, 0,3, D). Despite their objectivity, these coefficients need not reduce to scalars. 
In this case, the full system of governing equations is 

g = —gdivv, 

gc = —div3 + m, 
(jib 

gv = -grad(g 2 -^-{g,c)) +div(U(g,c,3,D)3 + V{g,c,3,D)D) +b, 

Z(Q,c,3,D)3=-grad(^(g,c))-Y(g,c,3,D)D+£, 

(3.19) 
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These equations are substantially more complicated than those for the case of simple 
diffusion without convection. However, if we require that h and T v - ls be linear functions 
of (j,D), the general result (3.17) specializes to yield 



h = -C(g,c)j, 

vis 



with 



C(e,c)>0, 7]x(g,c)>0, and 2r)i(g, c) + 3r/ 2 (g, c) > 
for all (g,c). In this case, the full system of governing equations has the form 

-gdivv, 



(3.20) 



(3.21) 



Q 

QC 

gv 

({S, c)3 

which constitutes a system of coupled equations for the fields (g,c,v,j 



-div j + m, 

-grad (g 2 -^(g, c)) + div (2ij 1 (g, c)D) + grad (772(0. c)divv) 
dip 

-grad( — (g,c)) +£, 
oc 



> (3.22) 



b. 



3.2.2 Incompressible mixture 

If the convective mixture velocity obeys 

divv — tr D = 0, 



(3.23) 



so that the motion of the mixture is isochoric, we may without loss of generality assume 
that g = 1 and consequently delete the mixture density from the list of independent 
constitutive variables. For convenience, we introduce the fields 

p mixture pressure, 
S mixture extra stress, 

defined by the decomposition 

T=-pl+S, p=-\trT, trS = Q, (3.24) 

of the mixture stress into a constitutively indeterminate pressure p that reacts to the 
constraint ( 3.23| ) and a traceless but constitutively determinate extra stress S. 
We then assume that 



with 



W>, /i, n, h, S) - (i>(z), A(«), II{z), h{z), S{z)), 



z = (c,j,D), 



(3.25) 



(3.26) 



and, requiring that these relations satisfy the energy imbalance ( |3.9| ) in all processes, 
find that the constitutive response functions ip, /t, and IT delivering the mixture energy 
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density -0, the solute diffusion potential fx, and the solute microstress 77 must, as in 
( 2.13 ), be independent of the diffusive mass flux j and the strain-rate D, so that 



and 



n 



4'{c)i. 



(3.27) 



Further, granted that the constitutive response functions h and S for the internal solute 
body microforce h and the traceless component S of the mixture stress are smooth 
functions of (c, j, D), we find that 



(3.28) 



h = -Z(c,3,D) 3 -Y(c,3,D)D,} 

S= U(c, 3 ,D)3 + V(c,3,D)D,j 
where Z, Y, U, and V must obey 

3 -Z(c, 3, D)j + 3-Y(c,3, D)D + D-U(c l3 , D)j + D-V{c, 3l D)D > (3.29) 
for all (c,j, D). 

If, as in the compressible case considered above, we require that h and S be linear 
functions of (j,-D), the general result ( 3.28| ) reduces to 

h = -C(c) 3 , 1 
S= 2 m (c)D,\ 

with 

C(c) > and r?i(c) > 

for all c. In this case, the full system of governing equations is: 

= divv, 
c = —divj + 771, 

■it = —gmdp + div (2rji(c)D^j + b, 
C(c)j = -grad (4>'(c))+£, 

which constitutes a system of coupled equations for the fields (p, c, v, 3). 



(3.30) 



(3.31) 



(3.32) 



4 Discussion 

Our theory is based upon the interpretation of the diffusive flux as a generalized velocity 
and the introduction of corresponding power- conjugate forces, which we call microforces, 
into an otherwise conventional continuum-mechanical description. A systematic, ther- 
modynamically consistent derivation generates a final system of governing equations — 



( 2.22 ) when convection is ignored, and ( 3.22 ) or ( |3.32| ) when convection is taken into 



account — that constitutes a generalization of common convection-diffusion equations, a 
generalization that accounts for constitutive nonlinearities and external forces acting on 
the diffusing constituents. 

In the absence of mixture convection, we find that the constitutive response of the 
material is determined by the provision of two scalar functions, one for the energy density 
and the other for the frictional resistance. This result generalizes naturally to account for 
dissipation associated with the convective motion of the mixture. Importantly, we arrive 



at our theory without introducing a constitutive relation such as (1.3) for the diffusive 
flux. 



16 



E. Fried & S. Sellers 



Additionally, our proposed format also leads to a generalization, ( [2.25 ), of the classical 



relation between diffus ion coefficient and frictional resistance. We find that the classical 
form for the relation (|l.9| ) holds when the free-energy density is given by the special 



constitutive relation (2.28) 



A distinction between the final governing equations of our theory and those arising 
from more standard derivations of diffusion equations is that, in our approach, a vector- 
valued microforce balance provides a generalization between diffusive solute flux and 
gradient of the solute diffusion potential. Furthermore, consistent with Nernst's treat- 
ment of osmotic forces, the gradient of solute diffusion potential emerges naturally as the 
diffusive force in this balance. This result for the diffusive driving force is consistent with 
Batchelor's |34|] "thermodynamic force" on particles in suspension. Typically, however, 
this diffusive force is not formally included in force balances, but, rather, a diffusive flux 
is added ad hoc to the total flux. We point out that Bird, Armstrong & Hassager (30) 
include such a diffusive force in their force balance for polymer suspensions. 
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